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Abstract It has been conjectured that the relative phase between strong and electromagnetic amplitudes is 
universally —90° in charmonium decays, ip' decaying into pseudoscalar pair provides a possibility to test this 
conjecture. However, the experimentally observed cross section for such a process is depicted by the two-fold 
integral which takes into account the initial state radiative (ISR) correction and energy spread effect. Using 
the generalized linear regression approach, a complex energy-dependent factor is approximated by a linear 
function of energy. Taking advantage of this simplification, the integration of ISR correction can be performed 
and an analytical expression with accuracy at the level of 1% is obtained. Then, the original two-fold integral 
is simplified into a one-fold integral, which reduces the total computing time by two orders of magnitude. Such 
a simplified expression for the observed cross section usually plays an indispensable role in the optimization of 
scan data taking, the determination of systematic uncertainty, and the analysis of data correlation. 
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1 Introduction 



The relative phase between the strong and the 
electromagnetic amplitudes of the charmonium de- 
cays is a basic parameter in understanding decay 
dynamics. Studies have been carried out for man; 
J/tp two-body decay modes: 1~0~ [l|, 3) 0~0 
1, 1-1- i and^fl. These analyses Reveal that 
there exists a relative orthogonal phase between the 
strong and the electromagnetic amplitudes in J/tjj 
decays [l|-0|- As to ip', there is also a theoretical 
argument which favors the ±90° phase 0]. Experi- 
mentally, some analyses based on limited 1 
and 0"0~ data indicate that the large phase is com- 
patible with the data. Moreover, some efforts have 
been made to extend the phase study to ip" decay 
phenomenologically 12, 3| and experimentally 



14 



The great merit of the phase study lies in that 
it can provide the valuable clue for the relation be- 
tween the strong and the electromagnetic interac- 
tions. Now with the upgraded accelerator and de- 
tector, BEPCII/BESIII, on May 2009, the high lumi- 



nosity of 3 x 10 32 cm" 2 s _1 had achieved, which is the 
highest luminosity in r-charm energy region which 
ever existed. The 106 M ijj' and 226 M J /tp events 
have been collected, even more colossal data are to 
be collected in the forthcoming years, which gives a 
great opportunity to determine the phase between the 
strong and the electromagnetic amplitudes with un- 
precedented statistical precision. 

A favorable way to measure the phase is through 
the scan experiment which is the most model- 
independent approach. However, even with high lu- 
minosity accelerator, the exclusive scan experiment of 
charmonium decay is fairly difficult due to low statis- 
tics at each energy point. Therefore, the optimization 
study for the data taking strategy is of great impor- 
tance in order to obtain the most accurate results 
with the limited luminosity (equivalently within the 
limited data taking time). 

Without losing generality, we focus on the mode 
of ijj' decays to two pseudoscalars. Because, as will 
be shown in the next section, this decay mode can ac- 
commodate a comparatively simple parametrization 
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form which is of great benefit to extract the relative 
phase. To get the optimized data taking scheme, we 
resort to the sampling simulation technique which is 
successfully used in the study of the data taking strat- 
egy for a high precision r mass measurement 
For such kind of method, great many times of fits 
should be done, where a large number of calculations 
need to be performed to get the theoretically expected 
observed cross section. Unfortunately, two nested in- 
tegrations in this calculation take too long time to 
make the actual optimization procedure impractical. 

In this paper, we devote to the simplification of 
calculation for the observed cross section of ip' decay- 
ing to pseudoscalar pair. Some reasonable assump- 
tions lead us to obtain the analytic expression for the 
Initial State Radiative (ISR) corrected cross section. 
That is to say, we transform the two-fold integral into 
a one-fold integral which speeds up the calculation by 
one hundred times. 

2 Observed cross section 
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where the real parameters <f> and C are the relative 
phase and the relative strength between the strong 
and the electromagnetic amplitudes, and B(s) is de- 
fined as @ 



B(s) : 
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Here y/s is the center of mass energy, a is the QED 
fine structure constant; M^,, and T t are the mass and 
the total width of ip'; T ee is the partial width to e + e~. 

The Born order cross sections for the three chan- 
nels read 



= s 3 / 2 
x \ J- „+. 
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The process of ip' decays to Pseudo-scalar and 
Pseudo-scalar (PP) final state can be parameterized 
by merely two amplitudes 0, E3] , that is 
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where A EM denotes the electromagnetic amplitude 
and A s the SU(3) breaking strong amplitude. Here, 
the G-parity violating channel tt + it~ is through the 
electromagnetic process (the contribution from the 
isospin-yiolating part of QCD is expected to be 
small [3] and is neglected), K° S K L through the SU(3) 
breaking strong process, and K + K~ through both. 
For e + e~ experiment, the actual amplitudes must in- 
clude the contribution of continuum which features 
the electromagnetic process 



A Tl 
A 



= A 



EM 



K+K- 
A K »K° 



= A 



Aem 
em -r Aem - 
As , 



-Ac 



(2) 



<(«) = 



^C 2 |S(.s)| 2 |^o K o( S )| 2 P K o K o( S ) 



(7) 

where C = 1+Ce 11 *; J"/. s .(s) = ff.s./s, with f f . s . 
being an energy independent constant, and f.s. — 
tt+tt-^K+K^K^KI; V f . s .(a) = 2q), J3s, with qj, 3 , = 

It is obvious that in Eq. © , if C = 1 , (t^+k- 0) is 
identical to cr^+^L (s) while if 0$ =Ce i *, o^+k- ( s ) is 
identical to cr^ r ™ (s). From the mathematical point 
of view, the cross section expression of &k°+k- ( s ) ^ s 
more general with the expressions of <J^+ r ™- (a) and 
a?, o r Il (s) as its special cases. Therefore, in the fol- 
lowing study, only tackled is the formula for K + K~ 
final state and f.s. is simply donated as K. 

In e + e~ collision, the Born order cross section is 
modified by the ISR in the way 
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where A C EM is the amplitude of the continuum con- 
tribution. Besides the common part, A C EM , A EM and 
A s can be expressed explicitly as 



where X f = l—s'/s. F(x,s) has been calculated to an 
accuracy of 0.1% 2^- 22 1 and H(s) is the vacuum po- 
larization factor. In the upper limit of the integration, 



2 



V? is the experimentally required minimum invariant 
mass of the final particles. In this work, X f = 0.15 
is used which corresponds to invariant mass cut of 
3.4 GeV/c 2 . 

By convention, r ee has the QED vacuum polar- 
ization in its definition 23J, l24[ . Here it is natural to 
extend this convention to the partial widths of other 
pure electromagnetic decays, that is 
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where 
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with vacuum polarization effect included. 

The e + e~ colliders have finite energy resolution 
which is much wider than the intrinsic width of ip'. 
Such e nerg y resolution is usually a Gaussian distri- 
bution 
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where W = y/s and A, a function of the energy, is 
the standard deviation of the Gaussian distribution. 
The experimentally observed cross section is the ra- 
diative corrected cross section folded with the energy 
resolution function 
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Based on Eq. ^ , the whole expression of the ob- 
served cross section is subdivided into three terms: 
the continuum, the resonance, and the interference 
terms. The simplification of each term will be dis- 
cussed separately. 

3.1 Continuum term 

In the light of Eq. ([6]) , the Born order expression 
for the continuum is written explicitly as 



a 

° Born 



8ira 2 f 2 K (s/A-m 2 K ) 3 / 2 



In the above equation, 
factor 
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the most crucial part is the 



'9/2(5) : 
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(T obs {W) = J dW'a r . c .(W')G(W',W). (10) 


For briefness, the variables r ee , M^,/, and T t are 
respectively written as r ee , M, and T hereafter. 

3 Simplification of ISR correction 

In this section, we focus on the simplification 
of ISR correction of the observed cross section. In 
the energy region we concerned (3.67 GeV/c 2 - 3.71 
GeV/c 2 ), the vacuum polarization factor could be 
concerned as constant and absorbed into T ee as in 
Eq. (|9]). So we could begin with this expression: 



0V.c(s) : 



dxF(x,s)a Born (s(l-x)) . (11) 



F(x,s) is the structure function, which can be ex- 
pressed as follows: 

F(x, s) = X 1 ' 1 ■ B^+x* ■ B 2 (t) + 

x t+1 -B 3 {t) + 0(x t+1 t 2 ) , (12) 
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Fig. 1. Variations of factor lg(s) against center- 
of-mass energy (\/s) in the vicinity of ip' res- 
onance peak for ,8 = 9/2, 4, 7/2. 

For the study of charmonium physics, s is much 
greater than m 2 K , therefore the factor £g/2(s) variates 
almost linearly in the vicinity of ip' peak as shown in 
Fig.[TJ With this observation, it is natural to approx- 
imate the factor Zg/2(s) with a linear function, viz. 



; 9/2 (s)~Ag /2 -S + C 9/2 
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As a matter of fact, the similar factors appear in 
the resonance and interference terms as well. So gen- 
erally, we define 



his)'- 



(S/4- 77^)3/2 



and utilizing the approximation 
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Here the coefficients A^ and can be determined 
analytically, the details are degraded into the ap- 
pendix 0. 

With the linearization of the factor l 9 / 2 (s), the x- 
concerned ISR integral for the continuum term has 
actually the form 



Po- 



x^dx , 
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which can be integrated easily. So the ISR corrected 
cross section of the continuum is expressed analyti- 
cally as follows 



with 
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3.2 Resonance term 



In the light of Eq. ©, the Born order expression 
for the resonance is written explicitly as 
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where 



A 1 = 9T 2 e Ja 2 -{l+C 2 + 2Ccosc 
As far as the factor 
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is concerned, the similar approximation as the previ- 
ous section is adopted, viz. 

^7/2(s) ~ A 7/2 • S + C7/2 • 



The x-concerned ISR integral for the resonance 
term then reads 
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which can be integrated analytically 2JJ, [28 1 
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With the expression of p(s,t), the ISR corrected 
cross section of the resonance is re-casted as 



lna 2 f 2 K 



A r [(A r / 2 -s+C7/2>G (s)-A r/2 -s-G 1 (s)] , 
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3.3 Interference term 



The Born order expression for the interference can 
be acquired readily from Eq. ([6]). However, for clear- 
ness the expression of the interference is further di- 
vided into two sub-terms as follows 

h _8ira 2 f 2 K A 2 -(s-M 2 ) ( s /4-m 2 ) 3 / 2 



and 
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A 2 = 6(r ee /a)-(l+Ccos^ 



A 3 = 6(T ee /a)-CMTmn(p . 



*The determination of linear coefficients \p and (^g is similar to that of linear regression, where the optimization is used. 
However, for linear regression, a linear function is used to fit a set of separated data while for our problem, a linear function is 
used to approximate another non-linear function. Such an idea of linearization is referred to as the generalized linear regression. 
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The simplification strategy is the same as those 
used for the continuum and resonance. First, the fac- 
tor 

(£/4-m|)^ 
'4 S J- ^ 

is approximated as 

i 4 (s)« A 4 -s+C4 ; 

second, the x-concerned ISR integrals for the interfer- 
ence terms have the forms as those in Eqs. (|17j) and 
(f20|) . which can be integrated out directly or by For- 
mula (f2"Tj) . Finally, the ISR corrected cross section of 
the interference is obtained 

= • A 2 • { (A 4 • ,s + C 4 )(s - ' Go 

- [2A 4 • s 2 + (C 4 - A 4 M 2 )s] • G_ (s) + A 4 s 2 • G 2 (s)} , 

(26) 

^ 2 c . = ^^^3-[(A 4S +C4)-Go( S )-A 4S -G 1 ( S )] , (27) 

where G M (s) is given by Formula (|23|) . 

In summary, the ISR corrected cross section for- 
mula is 

VrAs) = <ZAs) + °L(s) + tfM + <#M , (28) 

with expressions of the cross section for each term 
given in Eqs. (HU), ([22]), ([26]), and ([27]), respec- 
tively. 

4 Possible simplification for energy 
spread integral 

As indicated in Eq. (|10l) . the experimentally ob- 
served cross section is the cr r c . convoluted G(W, W), 
which might be simplified further. Two methods, the 
Taylor Expansion (TE) method and the Fast Fourier 
Transformation (FFT) method, have been considered 
for such a simplification. 

For the TE method, we begin from Eq. ([TO]), and 
Taylor expand the a r c . at W, viz. 

^ cr (n) (W) 

where a r n } (W) denotes the n-th derivative of func- 
tion cr r c . at value W. Replace the Taylor expansion 

^http: / /on. wikipedia.org/wiki/Fast_Fourier_transform 



of <7 r . c . into Eq. (ITU]) , the integral to be calculated has 
the following form 

J x"e~ x2 dx , 

— oo 

which can be precalculated. However, in order to 
achieve a reasonable precision, we need to calculate 
hundreds, or even thousands of terms in Taylor ex- 
pansion. This means that the fairly high order deriva- 
tives of oy. , have to be calculated, and too much time 
is consumed, which is not acceptable. 

As to FFT method.., we could easily find that 
the observed cross section <J obs (W) is a convolution 
of the radiative corrected cross section and a gauss 
function. Considering the Convolution Theorem in 
Fourier Transformation 

where $ represents Fourier Transformation, (^) rep- 
resents convolution. To calculate convolution effi- 
ciently, we use Fast Fourier Transformation. First, 
<7 r . c . and G should be sampled in energy region. After 
that, we get two series of numbers. Then DFT (Dis- 
crete Fourier Transformation) should be performed 
on both series, and the resulting series should be 
multiplied to generate one final series. Finally IDFT 
(Inverse Discrete Fourier Transformation) should be 
performed on this series and what we get is the distri- 
bution of a obs in energy region on which a r . c , and G 
are sampled. This process is very fast, and we could 
get the result on the whole energy region at the same 
time rather than calculating the integral one by one. 
To get accurate result, the sample number should be 
very large (512 or 1024), which means a large number 
of cross sections should be calculated. In real energy 
scan, the number of data taking points is usually not 
large (less than 20). The total integration time in a 
small number of energy points is less than the time 
cost by sampling a large number of cross sections and 
perform DFT and IDFT on it. So this method does 
not fit our purpose. 

5 Investigation of simplified formula 

5.1 Precision 

The accurate observe cross section (a obs ) is cal- 
culated by Eq. ([TO]) while the simplification one (de- 
noted by cr a obs ) is also calculated by Eq. ([TO]) but with 
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(ov. c .) replaced by the expression (f^5|) . The relative 
error of two observed cross sections is defined as 



Ra 



(29) 



In the calculation of the observed cross section, all 
parameters of resonances are taken from PDG08 [29j, 
A = 1.3 MeV is used. Two real undetermined pa- 
rameters are the relative phase (</>) and the relative 
strength (C) between the strong and the electromag- 
netic amplitudes. The dependences of R a on <p and C 
are shown in Figs. [5] and [3] respectively. 




3.67 3.68 3.69 3.7 ; 

Vs/GeV 

Fig. 2. Variations of R a against ^/s in the vicin- 
ity of i/i' resonance peak for = 0°, 90°, 180°, 
and 270°. In the calculation of the observed 
cross section, C is fixed at 2.5. 



c= 1 
c= 5 
c= 10 
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Fig. 3. Variations of R a against yfs in the vicin- 
ity of ip' resonance peak for C =1, 5, and 10. In 
the calculation of the observed cross section, 
<b is fixed at 90°. 



The variations of R a against the center-of-mass 
energy (y / s)in the vicinity of if>' resonance peak for 
= 0°, 90°, 180°, and 270° are displayed in Fig. [2j 



according to which we notice that firstly, the abso- 
lute value of R a is less than one percent in the en- 
ergy region we concerned; secondly, the difference be- 
tween two cross sections fades away at the resonance 
peak; thirdly, the differences in off-resonance region 
are larger than that in on-resonance region. The sim- 
ilar dependence of R a on C can be seen from Fig. [3l 
where displayed are the variations of R a against t/s 
in the vicinity of tjj' resonance peak for C =1, 5, and 
10. It is obvious that the difference due to the vari- 
ation of C is even smaller, which is at the level of a 
few per millc. 

5.2 Computation time 

The symbol T s (T°) denotes the computation time 
when a s obs (a obs ) is used for the cross section calcula- 
tion. The comparison of T s (denoted by the solid 
line) and T° (denoted by the dotted line) at both res- 
onance and off-resonance regions are shown in Fig. HJ 
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3.67 3.68 3.69 3.7 3.71 

Vs/GeV 

Fig. 4. Comparison of T s and T u at both reso- 
nance and off-resonance regions. 

From comparison it can be seen that about 
one-hundred-time reduction of computation time is 
achieved by our simplification algorithm. Although 
only one-fold integral is simplified by analytic expres- 
sion, the computation time is less than 0.1 second for 
each energy point which is fast enough for our scan 
simulation study. 

5.3 Application 

As we mentioned in the introduction, the speed of 
calculation of the observed cross section is the crucial 
issue of data taking optimization study of the scan 
experiment. Without reasonably simplified formula, 
it will be a too long time to perform the optimization 
fit, and the detailed scan optimization is impractical. 

Besides the application in scan optimization, sim- 
plified cross section formulas can also used for the un- 
certainty study [30( and correlation study [3l| . Since 
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for both of these studies, the sampling-and-fitting 
method is also adopted, the fast computation of cross 
section is needed as well. 

6 Summary 

The complete expressions for ip' — > PP decays are 
presented, including the relative phase between the 
strong and the electromagnetic amplitudes. After 
linearizing one non-linear kinematic factor, the in- 
tegrand with the initial state radiation is integrated 



analytically. Such a simplification of two-fold integral 
into a one-fold integral reduces the total computing 
time by about one hundred times. 

The possible approaches for simplification of en- 
ergy spread integral are also discussed. 

The simplified formulas of the observed cross sec- 
tions obtained in this paper provide a practical tool 
for the further optimization study of the scan data 
taking, which is of great importance for the study of 
the relative phase between the strong and the elec- 
tromagnetic amplitudes. 



Appendices A 

As we have noted in Subsection 3.1, the factor 



k/2(s) = ■ 



,9/2 



varied almost linearly in the vicinity of ip' peak, and its 
variation against s is shown in Fig. [T] Therefore, for the 
factor 



his) 



(s/4-m^ 
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a linear function (it refers to Eq. ()16p ) 



lp(s)^\ B -s + (p (A2) 

is utilized to approximate it in the vicinity of resonance 
peak. The coefficients \p and are determined by the 
generalized linear regression method. As the first step, we 
define the integration 



C 2 : 



Both Ci and C2 contain integral 

d(/W s V^ )3/2 

J S1 

where 
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/3 = 2,-,3,~,4,~ 
'222 

u = 4m 2 K . 

For different /3, we can calculate the integral analyti- 
cally 0. For P = 2, 



D(2) = 
For P=- 



\Jx — u ( — + 2 ) — 3\/wtan 
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The needed values of coefficients A^ and C,p are ob- 
tained by the minimization of the integration /, that is 



|L=0 and «0 



(A4) 



From the above requirements, we acquire a set of lin- 
ear equations of A^ and solve it, we obtain 



S1C1— S2C2 . . 83C2 — S2C1 
X f*= JTT. *2 and C/3 = f f r2 



where 
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For /3 = 3, 
D(3) = 
For /3 = 



3 _i / \fx — u 
tan 



4^ 



1 (2u 5 



4 \r 
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For p = 4, 



8m* 



For /? = - 
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2' 



Y M _L 7_ 

\3x 3 + 8ux 12x 2 

(All) 



-fThe following integrals are obtained by using Mathematica and checked by hands. 
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9\ 2(x-u)?(5u + 2x) 



D I - ) = ^_^r^^±l . (A12) 

It could be easily checked that for the coefficients \p 
and C,p we obtain, 



d 2 i 



—2=2/ ds = 2(s 2 -si)>0 . 



(AI3) 



(A14) 



This means what we get is the minimum of /, not maxi- 
mum. 

The relative error between the linearized formula and 
the original formula is defined as follows: 

\b-h\ 



Ri 



(AI5) 



When /3 = 9/2, 4, 7/2, the variations of 7?; against the 
center-of-mass energy (^/s) are shown in Fig. [5] 




Fig. 5. The variations of Ri against t/s for 
= 9/2, 4, 7/2. 



References 



1 Jousset J et al (DMII collaboration). Phys. Rev. D, 1990, 
41: 1389 

2 Coffman D et al (Mark III collaboration). Phys. Rev. D, 
1988, 38: 2695 

3 Suzuki M. Phys. Rev. D, 1999, 60: 051501 

4 Lopez G, Lucio M J L, Pesticau J. hep-ph/9902300 

5 Kopke L, Wermes N. Phys. Rep., 1989, 174: 67 

6 Baldini R et al. Phys. Lett. B, 1998, 444: 111 

7 Suzuki M. Phys. Rev. D, 2001, 63: 054021 

8 Gerard J M, Weyers J. Phys. Lett. B, 1999, 462: 324 

9 YUAN C Z, WANG P, MO X H. Phys. Lett. B, 2003, 567: 
73 

10 WANG P, YUAN C Z, MO X H. Phys. Rev. D, 2004, 69: 
057502 

11 BAI J Z et al (BES collaboration). Phys. Rev. Lett., 2004, 
91: 052001 

12 WANG P, YUAN C Z, MO X H. Phys. Lett. B, 2003, 574: 
41 

13 WANG P, MO X H, YUAN C Z. Int. J. Mod. Phys. A, 
2006, 21: 5163 

14 Ablikim M et al (BESII collaboration). Phys. Rev. D, 2004, 
70: 077101 

15 WANG Y K, MO X H, YUAN C Z et al. Nucl. Instrum. 
Methods A, 2007, 583: 479 



16 WANG Y K, ZHANG J Y, MO X H, YUAN C Z et al. 
China Physic C, 2009, 33: 501 

17 Haber H E and Perrier J. Phys. Rev. D, 1985, 32: 2961 

18 Chernyak V L and Zhitnitsky A R. Nucl. Phys. B, 1982, 
201: 492 

19 WANG P, MO X H, YUAN C Z. Phys. Lett. B, 2003, 557: 
192 

20 Kuraev E A, Fadin V S. Sov. J. Nucl. Phys., 1985, 41: 
466-472 

21 Altarclli G and Martinelli G. CERN, 1986, 86-02: 47 

22 Berends F A, Burgers G, Neerven W L. Nucl. Phys. B, 
1988, 297: 429; Berends F A, Burgers G, Neerven W L. 
Nucl. Phys. B, 1988, 304: 921 

23 Tsai Y S. SLAC-PUB-3129, 1983 

24 Alexander P et al. Nucl. Phys. B, 1989, 320: 45 

25 Lee S Y. Accelerator Physics (2nd Edition). ShangHai: Fu- 
Dan University Press, 2006 

26 Wille K. The Physics of Particle Accelerators. New York: 
Oxford University Press, 2000 

27 CHEN F Z,WANG P,WU J M, ZHU Y S. HEP & NP, 1990, 
14: 585-595(in Chinese) 

28 Cahn R N. Phys. Rev. D, 1987, 36: 2666 

29 Amslcr C et al. Phys. Lett. B, 2008, 667: 1 

30 MO X H, ZHU Y S. HEP & NP, 2001, 25: 1133-1139 (in 
Chinese) 

31 MO X H. ij}(2S) scan and some other studies (Post-doctor 
report). Beijing, 2003 



8 



